Thermomechanical squeezing below the 3dB limit 
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While cooling techniques are commonly used to remove thermal noise in an oscillator, parametric 
modulation of the oscillator can reduce both thermal and zero-point noise in one quadrature of 
motion. This "noise squeezing" property has been widely demonstrated, but has found little use in 
practice due to a well-known 3dB limit to the squeezing in equilibrium. Using a standard AFM can- 
tilever, we experimentally break this limit by combining optimal estimation and detuned parametric 
amplification. We observe squeezing of the noise by more than 5dB when the parametric modu- 
lation is applied. This squeezing — limited only by the oscillator Q-factor — could be used at room 
temperature to improve the stability and resolution of sensors, as well as in the low temperature 
regime to achieve mechanical squeezing below the quantum zero-point motion. 



High-quality mechanical oscillators are widely used for 
weak force detection[l|, Q, nano-scale controljj, Q and 
quantum state engineering [B], Periodic modulation of 
the spring constant or effective mass of a harmonic oscil- 
lator at twice its resonance frequency fo leads to para- 
metric amplification [7j, which enables a range of practical 
applications such as sensing[8j and information storage 
and processing 9] . In sensing, for example, mechanical 
oscillations that are in-phase with the drive can be pre- 
amplified above the measurement noise floor. This is 
commonly employed in nano-electromechanical systems 
where measurement noise is generally dominant 0, lioj . 
Such phase-sensitive amplification is also noiseless, which 
is useful in a quantum measurement context to observe 
signals below the standard quantum limit (SQL) on 
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These pre-amplified signals can then be 
transduced without requiring the strong measurement as- 
sociated with other means to surpass the SQL, such as 
back-action evasion[12|. 

The quadrature of motion out-of-phase with the para- 
metric drive exhibits a reduction in noise to below the 
thermal level 1(1 12 1 , lea ding to additional applications 
of parametric modulation |14l4l6| . This "squeezed" ther- 
mal distribution can, for instance, be used to mitigate 
the effect of other nonlinearities in a driven oscillator by 
reducing either phase or amplitude noise [l4j. Further- 
more, the noise in one quadrature can be further reduced 
to below the ground state variance V g , an effect referred 
to here as "quantum squeezing" . This kind of nonclassi- 
cal mechanical state is a precursor to tests of macro-scale 
entanglement and quantum gravity fl7| . However, in all 
previous experiments the steady-state squeezed variance 
has been constrained to a minimum of 50% of the orig- 
inal variance by a well-known 3dB limit [18]. Since the 
initial oscillator variance is Vr — V g (2N t h + 1), where 
Nth is the mean phonon occupation from coupling to a 
thermal bath, this limit leads to a stringent tempera- 
ture requirement of Nth < 1/2 f° r steady-state quan- 
tum squeezing. For typical micro- and nano-mechanical 
oscillators with resonance frequencies in the range of 1- 



FIG. 1. Schematic of the experimental setup. The red path 
denotes the high-fidelity carrier signal while the green path 
denotes the low-fidelity signal created by amplitude modula- 
tion (AM). 



100MHz, this imposes a pre-cooling temperature require- 
ment between 0.05 and 5mK, outside the range of most 
conventional cryogenic setups. Here, we demonstrate 
squeezing below 3dB for the first time, both transient 
and in the steady-state. The former is achieved by oper- 
ating outside of equilibrium, as demonstrated previously 
in ion traps [3], while conditional steady-state squeez- 
ing is demonstrated by using detuned parametric driving 
and optimal estimation [20j. Our results provide a path 
towards enhanced force sensitivity 2l| as well as precise 
control and robust quantum squeezing of mechanical os- 
cillators at high temperature in the absence of strong 
measurement [20] . 

The experimental setup is shown in Fig. [TJ The po- 
sition of a standard AFM cantilever is monitored us- 
ing a Mach-Zender interferometer in a balanced homo- 
dyne configuration, with a fiber tip used to focus the 
optical field onto the cantilever. The cantilever reso- 
nance frequency /o is modified from 9.6kHz to 12.5kHz 
upon applying a 450V potential difference between the 
cantilever and electrode; an effect due to the nonlin- 
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FIG. 2. Evolution of the squeezed and antisqueezed quadratures with a drive applied at t = 0. (a) Variance vs time traces 
for an on-resonance drive (green) and a below-threshold detuned drive (purple) for % = 22.5Hz with the cantilever tuned to 
fo — 12.5kHz. Solid lines are theoretical fits, while points show experimental statistics generated from 200 iterations of the 
drive turn-on. At each point in time, the quadratures are rotated so that the covariance (XY) — (X)(Y) fa over all iterations. 
A ring-up time of 2.5ms is chosen for the parametric drive to minimise impulse forces on the cantilever. Saturation due to the 
limit of fringe visibility can be seen for amplification above threshold, coinciding with a degradation of measured squeezing, (b) 
Theoretical (left) and experimental (right) variances as a function of detuning and time. Blue areas indicate squeezing below 
3dB. 



ear position dependence of capacitive energy in this 
geometry Q. This large DC offset increases the peak- 
to-peak frequency modulation \ f° r a given superim- 
posed AC voltage Q . A lock-in amplifier with a band- 
width much wider than the oscillator decay rate 7 al- 
lows the dynamics to be observed in quadrature phase 
space, which shows the oscillation amplitude and phase 
with respect to a reference frequency. The orthogonal 
quadratures X and Y then appear as components of the 
positions = ATcos(27r/dt) + y sin(27r/di) where we define 
an arbitrary reference frequency fd = fo + A. The para- 
metric effect from an alternating voltage at frequency 
2fd applied between the electrode and cantilever leads to 
a preferred phase of mechanical oscillation and hence a 
squeezed thermal distribution in the X-Y plane. 

When driven on resonance (A = 0), and with a 
strength above the instability threshold (x > 7), squeez- 
ing below 3dB can be achieved while the orthogonal 
quadrature is amplified indefinitely. Figure [5^ shows 
the time-evolution of the maximally squeezed and anti- 
squeezed quadrature variances measured in this regime, 
with x — 22.5Hz and 7 = 2Hz. Good agreement with 
theory is observed at short times (<20ms) with expo- 
nential growth of the amplified quadrature and classical 
squeezing approaching lldB in the orthogonal quadra- 
ture. However, the amplified quadrature saturates af- 
ter approximately 35ms, where the amplitude approaches 
the optical quarter-wavelength of 195nm and oscillations 
are no longer confined to the linear portion of the interfer- 
ence fringe. Crucially, a side-effect of this measurement 
nonlinearity is a severe degradation in observed squeez- 
ing well before saturation is apparent. Such limits to 



dynamic range therefore preclude the generation of all 
but transient squeezing above threshold. Nonetheless 
the strong squeezing observed reproduces the nonequi- 
librium squeezing observed in trapped ions[19j for the 
first time in a micromechanical oscillator, albeit in the 
classical regime. 

By detuning the parametric drive off resonance, the 
oscillator phase undergoes a net rotation with respect to 
the amplification axis, increasing the instability thresh- 
old to xth — \/ A 2 + "Y 2 [20I] . Consequently for the drive 
strength used here (\ = 22.5Hz), the oscillator is un- 
stable for detunings below A = 22.4Hz. When the de- 
tuning is increased further so that A > x, the phase- 
space trajectories form stable elliptical orbits, as illus- 
trated in Fig. The variances initially mirror this os- 
cillatory behaviour before relaxing to steady-state values 
in the long time limit, in quantitative agreement with 
theoretical modelling jl8|. The effect of increasing detun- 
ing on the transient statistics can be seen in Fig. [2Jd, 
where a dramatic change from monotonic behaviour to 
clear oscillations in the variance occurs at the thresh- 
old detuning. Notably, instantaneous squeezing below 
3dB is still possible below threshold, owing to the rapid 
drive turn-on. The final steady-state variances can be 
expressed with respect to the thermal variance Vt as 
V±/Vt = (It x/Xth)^ 1 - The squeezing limit of Vr/2 
is therefore — while not a constraint outside of equilib- 
rium — a fundamental one when in the steady-state. 

From the above observations, the benefit of a paramet- 
rically driven system in equilibrium would appear to be 
limited to enhanced readout in one quadrature and 3dB 
reduced variance in the other. However, these phenom- 
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FIG. 3. Lowering conditional variance via estimation, (a) Typical phase-space trajectory over a short time with a detuned 
parametric drive applied to the oscillator. The estimates {X ea t, Y es t} at a given time are calculated by filtering the low- fidelity 
data, and agree with a high-fidelity measurement (blue curve) to within an uncertainty (yellow ellipse). A feedback force Ffi, 
can be modelled by subtracting the estimates from the high-fidelity data, (b) Quadrature phase-space trajectories for 22.5 
second samples obtained from high-fidelity measurement (blue) and the residual after subtracting the estimate (yellow). The 
upper panels show the random-walk pattern in the undriven case mixed down at the resonance frequency fo = 14.5kHz for 
weak (A), intermediate (B) and strong measurement (C). The lower panels show the elliptical trajectories and residual noise 
for a parametric drive strength of \ = 57Hz detuned close to threshold and using the same SNR as above. 



ena can be combined using detuning for precise position 
estimation [18|. For an oscillator that is detuned so that 
A > x and has relaxed to the steady-state, the thermally 
excited oscillations will alternate between amplified and 
squeezed quadratures before decaying. Since the dynam- 
ics of the system are well known, a measurement of the 
amplified quadrature will provide some capacity to es- 
timate the squeezed quadrature at a later time. The 
squeezed quadrature therefore obtains an effective sen- 
sitivity enhancement without amplifying its mechanical 
fluctuations. 

Using control theory, the filter that extracts the best 
quadrature estimates from noisy data can be writ- 
ten in terms of parameters defining the oscillator's 
expected time-evolution along with the measurement 
sensitivity [Tij. The effective precision of these estimates 
can then be quantified using conditional variances. For 
example, the conditional variance of the X quadrature 
is defined as Vx = ((X(t) - X est (t)) 2 ) where X est (t) is 
the best estimate available from all prior measurement 
results at time t. This quantity equals the unconditional 
variance {X 2 ) in the weak measurement limit, and is 
quantum-limited to the mechanical ground-state variance 
for an ideal continuous measurement without paramet- 
ric driving [ll|. The conditional variance also defines the 
minimum effective temperature of a quadrature achiev- 
able from ideal feedback cooling, equivalent to applying 
phase-space displacements in the X and Y quadratures 
by X est and Y est respectively over time. 

In order to verify that the expected filter parameters 



produce the best possible estimate and quantify the re- 
sulting conditional variance, we use a bright optical field 
to provide a high-fidelity measurement of X and Y. As 
shown in Fig. [TJ a weak sideband is created using an in- 
tensity modulation of the bright field in order to provide 
low-fidelity measurements X and Y with independent 
shot-noise characteristics. The estimates are obtained 
from X and Y in post-processing and compared to the 
high-fidelity measurements X and Y. This use of two 
independent measurements to calculate the conditional 
variance is depicted in Fig. [3p,. 

The strength of the low-fidelity measurement can be 
naturally expressed in terms of a signal-to-noise ratio 
(SNR), quantified here as the ratio of mean-square ther- 
mal displacement to the square of the minimum distance 
resolvable over a time At w I/47 22]. The SNR was set 
by varying the modulation depth and thus the sideband 
intensity. The two continuous measurements of the can- 
tilever motion were then recorded in both the undriven 
case and using a parametric drive \ — 57Hz, with a de- 
tuning of A = 63Hz to ensure the below-threshold condi- 
tion. Optimal X and Y estimates were then generated in 
post-processing by minimising the respective conditional 
variances over the filter parameters. 

For an undriven oscillator, the optimal filter is of the 
form X est (t) = g X(t) * e- Tot and Y est {t) = g Y(t) * 
e -r t j n casej the conditional variances should be 
equal, i.e. Vx = Vy = Vo and their reduction from the 
thermal variance Vr is a function of the measurement 
strength. Phase-space Brownian trajectories {X, Y}, de- 
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FIG. 4. Steady-state conditional squeezing. Squeezing ratios are plotted against SNR for A = 63Hz (left) and against detuning 
for SNR ~ 1 (right). Blue points show the conditional squeezing ratio, with theoretical fits shown as solid lines. Dotted red 
curves are fits to the unconditional squeezing, limited to 1/2, with shaded bands to represent the experimental error margin. 
Labels A-C indicate the datapoints used to generate the trajectories in Figure [3] 



termined from the high-fidelity measurement, and corre- 
sponding residual noise {X — X est , Y — Y est } after apply- 
ing this optimised filter are plotted in Fig. for weak, 
intermediate and strong measurements. As the thermal 
signal increases towards the noise floor (SNR approaches 
1) the optimal filter bandwidth T = 7(1 + 4SNR) be- 
gins to widen and the conditional variances drop sharply 
as expected. At maximum sideband intensity, the RMS 
uncertainty is reduced from the thermal value of 240pm 
to a conditional value of 60pm. 

With a detuned parametric drive turned on, the two 
quadratures are no longer independent and the estimates 
require a more complex convolution 
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where for x < A, the optimal filter matrix H(i) takes 
the general form[18| 
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Here g n , Q, and <fi are positive real numbers. These 
parameters, as well as the angle a of the maxi- 
mally squeezed quadrature, are functions of SNR, drive 
strength and detuning. As in the undriven case, the 
squeezed variance V a is obtained by optimising over all 
filter parameters and the squeezing angle a. 

The squeezing ratio V a /Vo determined from this anal- 
ysis is shown in Fig. 2] as a function of SNR and agrees 
well with theory. The obtained filter parameters are close 
to those derived from theory 18[ , as shown in the Supple- 
mentary Information. As expected, the variances repro- 
duce the unconditional squeezing in the weak measure- 
ment limit and the parametric drive has no effect in the 
strong measurement limit. However, in the intermediate 
regime where SNR 1 there is a distinct minimum of 
V a /Vo ss 0.3, breaking the 3dB limit by a significant fac- 
tor. As can be seen in Fig. @] (right), the squeezing can be 



improved further by adjusting detuning closer to thresh- 
old. These results can be understood by the fact that the 
effective increased sensitivity due to the parametric drive 
is of greatest benefit near the noise floor and with max- 
imal amplification of the orthogonal quadrature. Since 
the maximum squeezing is proportional to •y/x/T) it can 
be enhanced by increasing the parametric drive strength. 
In principle, this squeezing can be made arbitrarily large, 
subject to the condition x ^ fa- Since the required filter 
(Eq. [2J can be integrated into a feedback loop, this al- 
lows arbitrary suppression of one quadrature of motion, 
circumventing the usual constraint associated with the 
measurement precision. 

Although demonstrated with thermal fluctuations, our 
technique applies in the same manner to the zero- 
point motion of an oscillator, with the maximum re- 
duction in conditional variance V a /Vo independent of 
temperaturefi"i|| . The action of the quantum modifica- 
tion at low temperatures — known as backaction noise — 
is instead to limit Vq to be no lower than the ground 
state variance. Therefore our approach could enable 
strong conditional quantum squeezing and, with the ap- 
plication of feedback, the same level of unconditional 
squeezing}^. Experiments with sensitivity near the 
standard quantum limit (where SNR « 1 at zero tem- 
perature) have been recently performed in mechanical 
oscillators (25|. While purely measurement-based 
schemes exist to create mechanical squeezed states 12j, 
Hfjj], significant squeezing requires high measurement 
strength and efficiency such that SNR > (2N t h + l) 2 . 
This regime is yet to be demonstrated in mechanical os- 
cillators. Nanoelectromechanical systems are, however, 
commonly integrated with a parametric drive and can 
be pre-cooled to near the ground state[i3|. Such emerg- 
ing systems are therefore good candidates for quantum 
squeezing below the zero-point motion using our tech- 
nique, even without significant advances in transduction. 

We have observed parametric squeezing of a microme- 
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chanical oscillator exceeding 3dB for the first time, both 
transient and in equilibrium, with the latter breaking a 
well-known limit for parametrically driven systems. This 
demonstrates that the 3dB limit to steady-state para- 
metric squeezing is not fundamental, and facilitates the 
wider use of thermomechanical squeezing in numerous 
applications such as ultra-sensitive local force detection, 
stabilisation and precision control. The combination of 
parametric driving, measurement and estimation sheds 
light on the important interface between quantum mea- 
surement and control that is being approached most no- 
tably in opto- and electromechanical systems. The tech- 
niques introduced, if applied to state-of-the-art readout 
techniques and high quality oscillators, also open the 
door towards engineering of nonclassical states of meso- 
scopic mechanical systems. More broadly, these results 
demonstrate that combining nonlincarity with control 
can both overcome fundamental limitations on nonlinear 
processes, as well as localise mechanical motion beyond 
constraints imposed by the measurement sensitivity. 
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Supplementary Information: Filtering of time-series 
steady-state data 



o 10 
O) 

o 1 



O 0.1 
100 



N 

X 10 



o 

-° 















60- 

N 

X 40- 

a 



0.2- 



o 



UU — - I W a u w t 

20 ■ °° o* o 










o 








0.01 


0.1 1 10 

SNR 


100 



FIG. Al. Optimal filter parameters. Solid lines indicate the- 
oretical parameters for \ — 57Hz,7 = 2.6Hz and A = 63Hz. 
Circles indicate parameters obtained from post-processing. 



At each sideband intensity featured in Fig. |U simulta- 



neous 45-second traces of the lock- in outputs X,Y,X,Y 
were recorded using a data acquisition unit with no para- 
metric drive applied to the cantilever. Throughout the 
experiment, during these thermal measurements the os- 
cillator resonance fo was kept within 5Hz of the lock-in 
reference frequency. Data acquisition was then repeated 
using parametric drives of various detunings at the same 
sideband intensity. The lock-in reference frequency - 
now shifted by A — is kept phase-locked to the drive 
voltage. A lock-in time constant of r c = 300/xs was used 
so that the output oscillations (limited in frequency to 
A) are contained within the output bandwidth. 

In post-processing, the conditional variance Vx was 
found by convolving the time series with a filter function 
(truncated to 22.5 seconds) and calculating ((X — X est ) 2 ) 
over the second half of the data. In the thermal case, the 
filter is characterised by a gain go, a decay rate To and 
a small rotation frequency Oo to account for drifts in f Q . 
The optimum conditional variance was then found by 
minimising computationally over all filter parameters. 

In the parametrically driven case, the measurement 
quadratures were initially rotated to be aligned with the 
unconditional squeezing. Optimal estimates of the filter 
parameters were then calculated based on known quan- 
tities: thermal signal-to-shot noise level (SNR), decay 
rate 7, drive strength \ an d detuning A[18(]. These ini- 
tial parameters are shown as solid lines in Fig. IA1I for 
the dataset analysed in this paper. Since the conditional 
squeezing angle is a function of SNR and other parame- 
ters, a quadrature rotation a. is included among the pa- 
rameters. The obtained optimal filter agrees reasonably 
well with theoretical estimates, although more likely to 
converge to consistent values in the higher-fidelity regime 
and in regions with more distinct minima. This demon- 
strates that the optimal filter parameters described in 
Ref. [IH agree quantitatively with experiment. 



